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(57) Abstract 



A blood glucose monitonng systwn (100) includes a broadband light source (110) and a specially fabricated optical filter (120) for 
modulating opncal radiation to be transmitted through a fleshy medium. Optical radiation (125) which passes through the fleshy medium 
(130) is detected by an optical detector (140) which generates an electrical signal indicative of the intensity of die detected light Digital 
signal processmg is performed on the electrical signal to extract those optical characteristics of the electrical signal due to the optical 
characteristics of Ae filter ^d constiments of the fleshy medium other than blood glucose concentration. The monitoring system employs 
a unique double-log transfonnation to minimize errors due to indeteraiinaie path length variations of the optical radiation through the 
fleshy n^ium The monitoring system (100) further employs specialized signal processing to avoid inaccuracies due to the previously 
umdentified solvent effect which arises when glucose is dissolved into water. previouwy 
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BLOOD GLUCOSE MONITORING SYSTEM 
Background of the Invention 

Field of the Invention 

The present invention relates to noninvasiye systems for monitoring blood glucose and other blood 
5 constituent concentrations. 
Oesfffiption of th? BpjqfBfi Art 

In the past, many systems have been developed for monitoring blood characteristics. For example, devices 
have been developed which are capable of determining such blood characteristics as blood oxygenation, ghicose 
concentration, and other btood characteristics. However, significant difficulties have been encountered when 

10 attempting to determine blood glucose concentration accurately using noninvasive blood monitoring systems. 

The difficulty in determinmg blood ghicose concentration accurately may be attributed to several causes. 
Fvst, blood glucose is typically found in very low concentrations within the bloodstream (e.g., on the order of 100 
to 1,000 times lower than hemoglobin) so that such low concentrations are difficult to detect noninvashrely, and 
require a very high signal-to-noise ratio. Second, there has been a lack of recognition of the kinds of noise and the 

15 proper method to use when removing this noise. For example, noise can be classified as deterministic (definable) 
or stochastic (random) where either of these kinds of noise could be linear (added) or modulated (multiplied). 
Knowledge of the distinction between the various kinds of noise is essential for purposes of using the proper method 
of removing noise. Additionally, the optical characteristics of glucose are very similar to those of water which is 
found in a very high concentration within the blood. Thus, where optical monitoring systems are used, the optical 

20 characteristics of water tend to obscure the characteristics of optical signals due to low glucose concentration withm 
the bloodstream. Furthermore, since each individual has unique blood properties, each measurement typically requires 
cafibration for the particular individual. 

In an attempt to accurately measure blood glucose levels within the bloodstream, several methods have been* 
used. For example, one method Involves drawing blood from the patient and separating the glucose from the other 

25 constituents within the blood. Although highly accurate, this method requires drawing the patient's blood, which is 
less desirable than noninvasive techniques, especially for patients such as small children or anemic patients. 
Furthermore, when blood glucose monitoring is used to control the blood glucose levej, blood must be drawn three 
to six times per day, which may be both physically and psychologically traumatic for a patient. Other methods 
contemplate determining blood glucose concentration by means of urinalysis or some other method which involves 

30 pumping or diffusing blood fkitd from the body through vessel walls. However, such an analysis tends to be less 
accurate than a direct measurement of ghicose within the blood, since the urine, or other blood fhiid, has passed 
through the kidneys. This problem is especially pronounced in diabetics. Furthermore, acquiring urine samples is 
often inconvenient. 

Another proposed method of measuring blood glucose concentration is by means of optical spectroscopic 
35 measurement In such devices, light of multiple wavelengths may be used to illuminate a relathrely thin portion of 
tissue, such as a fingertip or an earbibe, so that a spectrum analysis can be performed to determine the properties 



I 



wo 96/41151 

PCT/US96/08S06 

•2- 

of the blood flowing within the iyuminated tissue. Although such a method is highly desirable due to itt nonin«ish« 
character and its convenience to the patient probtems are associated with such methods due to the difficulty m 
bolatbq, each of the etements within the tissue by means of spectroscopic analysis. The difficuhy in determtamg 
blood glucose concentration is further e»e«bated due to the low com:entration of ghicose within blood, and the fact 
5 that ghicose in bh>od ha. »ery similar optical characteristics to water. Thus, it is very difficult to distinguish the 
spectral characteristics of glucose where a high amount of water b also found, such as in human blood. 

As is wen known in the ert, different motecutes, typicaHy referred to as constituents, contamed within the 
medium have different optical characterbtics so that they are mora or less .bsort.ent et different wavelengths of 
m. Thus, by analyzing the characteristics of the fleshy medium at different wavelengths, an indicetion of the 
10 composition of the fleshy medium may be determined. 

Spectroscopic analysis is based in part upon the Beer-Lambert law of optical characteristics for different 
elements. Briefly. Beer-Larabert's law states that the optical intensity of Hght through any medium comprismg a 
smgte substance is proportional to the exponent of the path lengths through the medium times the concentration of 
the substance within the meiKum. That is, • 
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(1) 



where pi represents the path length through the medhim and c represents the concentration of the substance within 
tha medium. For optical media which have several constituent substances, the optical intensity of the light receded 
from the ilhmdnated medium wiH be proportional to the exponent of the path hngth through the medmrn time, the 
concentratron of the first substance tbnes an optical absorption coefficient associated ivith the first substance plus 
the path length times the concentration of the second substance times the opticel ebsorption coefficient associated 
with the second substance, etc. That is. 
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(2) 



where «C represents the opticel absorption coefficient. 

Summary of t he Inventlnn 

Due to the parameters required by the Beer-Lambert law. the difficulties in detecting glucose concentration 
ansa from the difficulty in determining the exact path length through a medium (resulting from transforming the multl- 
path »gnal to an equhralem single-path signal), as well as difficulties encountered due to lo^ signal strength resultant 
fr«n a low concentration of blood gtocose. Path length through a medh.m such as a fingertip or earlobe Is very 
diffttult to determine, since not only «e optical wavelengtits absorbed by the fteshy medium, but else the signals 
are scattered wnhin the medium end transmitted titrough different path iengti.s. Furthermore, a, mdhratiid by ti.e 
above equations. ^ measured signel IntiMisity does not vary linearly witi, respect to tite path lengtit Themf or. 
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variations in path length of multiple paths of fight through the medium will not result in a Gnear averaging of the 
multiple path lengths. Thus, it is often very difficult to determine an exact path length through a fingertip or earlobe. 
In addition to these difficulties, it has been found that there is significant difficulty in detecting glucose within water 
based upon Beer-Lambert's law. Specifically, It has been found that inaccurate measuremenu are often taken of 
glucose concentration within water when optical measuring instruments are calibrated so as to assume maxinuim 
transmission of optical wavelengths wiQ occur through pure water without glucose. The present inventors have found 
that glucose together with water absorbs less than pure water for certain absorption bands and absorbs more for 
other bands. 

In addition to the aforementioned difficulties, current optical spectroscopic devices, as identified by the 
inventors for use in the present invention, often require expensive custom-made filters which are used to generate 
a pattern of optical signals to be transmitted. One such filter, comrnoniy known as a dichroic filter, comprises a 
rotating optically coated disk which Includes regions of varying optical thickness. The regions on the dichroic filter 
are formed in a pattern so that rotation of the optical disk results in the transmission of selected optical bands. 
The high precision necessary for optically coating the filter substrate with various thicknesses of optical material on 
minute portions of the optical disk typically makes this coating process highly expensive. The present invention 
decreases the cost of a rotating dichroic filter by a factor of approximately 100 times by relaxing the specifications 
of the filter and compensating for the relaxation of filter specifications through more intensive signal processing 
steps. The filter constructed in accordance with the present invention allows from 10 to 100 times as much light 
to pass while maintaining the same resolution. 

One aspect of the present invention involves a system for non invasively monitoring blood glucose 
concentration within a patient's bloodstream. The system has a light source which emits optical radiation at a 
phirality of wavelengths. A receptacle recehres a fleshy medium of the patient and an optical detector is positioned 
to recehre light from the light source and attenuated by the fleshy medium. The optical detector is responshre to 
optical radiation of at least the plurality of wavelengths to generate an output signal indicative of the Intensity of 
the optical radiation. A signal processor is coupled to the detector to receive the output signal. The signal processor 
is responsive to the output signal to isolate portions of the output signal due to optical characteristics of the fleshy 
medium to provide a set of frequency response values. The signal processor has a linearization module which 
Gnearizes the set of frequency response values and analyzes the linearized data to determine the concentration of 
glucose within the patient's bloodstream. In n one embodiment, the linearization module comprises a double logarithm 
operation. 

In one embodiment, the light source comprises a plurality of emitters, each emitter transmitting light at a 
selected one of the plurality of wavelengths. In another embodiment, the light source comprises a broadband light 
source and the system further has an optical filter which selectively passes ones of plurality of wavelengths present 
in the broadband light source. 

In one embodiment, the detector comprises a single detector responsh/e to the ones of the plurality of 
wavelengths to provide and output signal indicathre of the sum of the intensities of the ones of the plurality of 
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wavatengths. tn another embodiment, the optical detector comprises a plurality of detectors each detector 
responsive to at least one of the plurality of wavelengths to generate an output signa|imdicathre of the intensity of 
the at least one wavelength. 

Another aspect of the present invention involves a non-invaslve blood qIucoso monitorinfl system whicb 
5 analyies blood glucose within a fleshy medium conteming blood. The system has e light source end an optical 
detector responsive to light from the Kght source to gemirate an output signal A compression device Is configured 
to provide for physical perturbation of the fleshy medmrn to «(pr.ss fluid from the f teshy medium. A signal processor 
» msponsne to a first output agnal from the optical detector when fluid is expressed from the f teshy medium «,d 
respons«e to a second output signal from the optical detector when fyd is not expressed from the fteshy medium 
10 to isolate information relating to the concentration of glucose in the blood. 

Yet another aspect of the present invention involves a method of non-invasively. determining blood ghicose 
concentration. The method involves a number of steps. A set of values indicative of optical characteri«ics of 
signrficant blood constituents is generated and a fleshy medium having blood is illuminated with light of a pluraOty 
of wavelengths. The ight is detected after attenuation of the light by the fleshy medium. A signal is generated 
from the detected Sght which is indicathre of optical characteristics of the fleshy medium. .Components of the signal 
whch are indicathre of the concentration of glucose in water within the blood are' isotated in response to the 
detected light and the set of vahies indicathre of optical characteristics of the significant blood constituenu. A vahie 
« then generated indicative of the gbcose concentration within the blood. In one embodiment, the sffilbM itad 
constituents comprise water, hemoglobin, oxyhemoglobin and glucose dissohred in fhnd. 

In one embodiment, the signal is iinearired to provide an indication of the concentration of blood ghieose 
Advantageously, the indication of the concentration does not vary significamly with respect to the path length 
through the medium. In one embodiment, the UnearizinB invohres generating a first set of values indicative of the 
optn:.! Characteristics of the medhm,, taking a first logarithm of the first set of values to generate a second set of 
vehie. and taking a second logarithm of the second set of vah.es to obtab. a linearized set of values that are 
«d«atwe of the concentration of the blood constituent. In one embodiment the Bnearization further involves 
transformmg tiie first set of vah.es usmg a polynomial equation to provide a transformed first set of vahies. 
These end other aspects are described herein in further detail. 

Brief Descri Dtion of th|; nrawip ^f 
Figure 1 is a schematic block diagram Ulustrating the main functional efaimi^s of the blood ghicose 
30 monitoring system of the present invention. 

Figure 2 is a flow dngram which very generally depicts the method employed in accordance with the 
present invention to obtain blood ghicose concentration. 

Figure 3 Is a flow diagram which ilhistr ates in greater detaa the method of th. present invention used to 
determine blood ghicose concentration within a patient. 
35 Figures 4A4C are graphs which depict the optical transmisshm characteristics of a dKhroic filter over 

different degrees of rotation and at three separate wavelengths. 



20 



25 



9 



WO 96/41 151 PCTAJS96/08506 

-5- 

Figure 4D Qlustrates a matrix used to specify the optical characteristics of the dichroic filter used in 
accordance with the present invention. 

Figure 5 depicts the optical transmlcivity of a fleshy medium, such as a finger, plotted against wavelength, 
the optical transmlcivity of a selected portion of the dichroic filter used in accordance with the present invention 
5 plotted against wavelength, end the combined optical transmichrity of the fflter and the patient's finger plotted 
egainst wavelength. 

Figure 6A illustrates a plot of optical transmission through a medium versus a product of the absorption 
coefficient associated with the medium, the concentration of a substance within the medium, and the path length 
through the medium. 

10 Figure 6B is a plot of a "near-log" of the transmission through an optical medium versus a product of the 

absorption coefficient associated with the medium, the concentration of the substance within the medium, and the 
path length through the medium. 

Figure 6C illustrates a plot of the second log of the transmichrity through an optical medium versus the 
logarithm of the product of the absorption coefficient of the medium, the concentration of the medium, and the path 
. 15 length through the medium. 

Figure 7 is a flowchart which depicts the method used in accordance with the present invention to perform 
mittalizatlon and calibration functions. 

Figure 8 is a flowchart which depicts in greater detail the method used in accordance with the present 
invention to linearize values within the finger optical response matrix as indicated within a subroutine block of Figure 
20 3. 

Figure 9 is a data flow diagram which schematically represents the processing steps performed on an 

ft *i 

optical s'^nai m order to determine glucose concentration based upon the illumination of a fleshy medium by the 
optical signal 

Figure 10 depicts the method employed in accordance with the present invention to manufacture the dichroic 
25 filter depicted in Figures 1 and 9. 

Figure 11 is a graph which plots glucose concentrations versus wavelength with average water 
concentration subtracted. 

Rgure 12 schematicafly depicts a hand*held glucose monitor such as may be constructed in accordance with 
the present tnventran. 

30 Figure 13 depicts a dichroic filter as constructed by conventional methods. 

Detailed Desc riotion of the Invention 
OVERALL BLOO D GLUCOSE MQMITORING SYSTEM 

Figure 1 is a schematic block diagram which depicts the main functional and structural elements of a blood 
glucose monitoring system 100 in accordance with the present invention. As depicted in Figure 1, the ghicose 
35 monitoring system 100 includes a broadband light source 110 for emitting optical radiatibh 115 over a broad band 
of optical wavelengths. In one advantageous embodiment, the optical source 110 comprises a 3,000-kelvm, quarU, 
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halogen lamp. In addition to the halogen lamp, the light source 110 may further comprise, in one embodiment, a 
first indium gallium arsenide (InGaAs) fight emitting diode (LED) which emits light at a wavelength of approximately 
1,300 nanometers, a second InGaAs LEO which emits fight at epproximately 1.060 nanometers, and a third InGaAs 
LEO which emite light at approximately 940 nanometers. 

In an altemathre embodiment several LED's having wavetengths tuned to stiiected frequency bands could 
be used in place of the lamp 110. Of coarse, such an embodiment is typically expensive due to the cost of multiple, 
specianyfabricated LED's. Tbus, the preferred embodiment includes e device for producing multiple waveteneths m 
a cost effective manner as described below. 

The optical radiation 115 emitted from the light source 110 is focused via a lens assembly 117 (which may 
comprise fiber optics or the fike) and passes through a filtering element 120. The filter 120 may. in one 
advantageous embodiment, comprise a diehroic fQter constructed in accordance with the teachings of the' present 
invention, as described below with reference to Figure 10. The diehroic fiher 120 comprises an optically transmissive 
rotatable disk substrate which is layered with optical coatings having different thicknesses so as to modulate the 
broadband optical radiation 115 through a spectrum from the near infrared (NIR) (e-g.; 700 nm) to the infra^d |IR] 
{e.g, 1.400 nm). The filter 120 further indudes an optically opaque strip 122 which may. for exemple. comprise 
brass or some other metal which is deposited radially outward from the center of the rilt« disk 120. The opaque 
strip provides a "0" location indicator. The fBter disk 120 is driven in a circular motion (as described more fully 
below) by a stepper motor in one preferred embodiment. Filtered opticel radiation 125 passes from the filter 120 
through a fleshy medium, perfused with blood such as a finger tip 130. In some applicstions, it may be desirable 
20 to provide a focusing lens between the filter 120 and the finger 130. 

A pressure application device 129 Is shown schematically in Figure 1 to surround the finger tip 130. As 
wlH be discussed in greeter detail below, the pressure appHcation device 129 is used to squeeze fluid out of the 
finger tip 130 for purposes of enhancing signal detection. 

Light 135, which passes through the finger 130, is mcident upon en optical detector 140. In one 
advantageous embodiment, the optical detector 140 comprises an InGeAs optical detector which has a wall-defmed 
optical response for wavelengths in the spectrum of interest in the preferred embodiment (e.g, between 8S0 and 
1,700 nm). : ' 

The output of the optical detector 140 connects to a pre-ampBfier 150 which' edvantageously comprises 
a class A finear amplifier having an upper dynamic Gmit of between 1 and 5 volts. The output of the pre-amp 150 
eomiecu to a bandpass filter 160. The bandpass filter 160 advantageously comprises a linear RC filter having a 
lower cut-off frequency of approximately 0.5 Hz and an upper cut^iff frequency of approximately 10 kHz. The lower 
and upper cutnjff frequencies have roil-offs of approximately 6 decibels per octave. 

The bandpass filter 160 comieeu to an analog-to-dlgital converter 170. which, in preferred embodiment 
comprises a A-Z converter that converts the analog electrical signal output from the bamlpass filter to a digital 
sflnaL The analog-to-digital converter 170. for example, may comprise e 16-bit. 20 kHz conversion rate 
analog-to-digital converter. One example is Model No. CS5317, available from Crystal Semiconductor. The digitized 
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signal which is output from the convenor 170 is then provided to the input of a digital signal processor 180 for 
signal processing. In one enibodiment the digital signal processor may be implemented in software within a 
computer. For example, an INTEL 486 MP, or an ANALOG DEVICES OSP chip. Model No. 21062, are two examples. 
The digital signal processor 180 outputs a value indicative of the blood glucose level of the blood within the fmger 
5 130. 

In operation, when light 115 is emitted from the broadband light source 110;lDVer a wavelength range of 
approximately 700 nanometers to 1,700 nanometers, this broadband light 115 shines through the rotating dichroic 
filter 120. It should be noted that the tight 115 is focused onto a portion of the filter 120 by means of fiber optics, 
a lens assembly (e.g., the lens 117), or the like. As the dichroic filter 120 rotates, the broadband light 115 is 

10 filtered through a portion of the dichroic filter 120 producing the filtered optical radiation 125. As indicated above, 
the dichroic fitter 120 is coated with optical layers of varying thickness so that different portions of the dichroic 
filter 120 pass different wavelengths of light. Thus, as the filter 120 rotates, the optical radiation 125 output from 
the filter includes optical radiation of various wavelengths. In one embodiment, a fiber optic is used to ooupla the 
optical radiation 125 emitted from a portion of the fater 120 to the petient's finger 120. It should be noted here, 

15 that since the optical characteristics of the filter 120 can be carefully measured and the rotational speed of the 
dichroic fiKer 120 is known, the pattern of optical radiation 125 emitted from the filter, 120 to illuminate the finger 
130 is well defined, and therefore, may be used during signal processing to determine the amount of attenuation 

i 

which is due to the optical filter 120. 

The optical radiation 125 which is used to illuminate the finger 130 passes through the finger 130 to 
20 produce the detectable light 135. As is wefl known in the art, some of the optical radiation 125 passes unimpeded 
through the finger 130, some of the optical radiation 125 is reflected within the finger 130 to produce scattering. 
The scattered radiation which is Uansmitted through the finger 130, together with the light which passes unimpeded 
through the finger 130, make up the light 135. Some of the optical radiation 125 is absorbed by constituents within 
the f inger 130. 

25 The finger 130 is known to include a fingernail, skin, bones, flesh, and blood. The blood itself primarily 

comprises water, oxyhemoglobin, reduced hemoglobin, lipids, protein and glucose. Each of these constituents within 
the finger (e.g., nerves, muscle tissue, etc.) contribute to the absorption and scattering of the optical radiation 125 
through the finger 130. As described above, the absorption of optical radiation through a nonhomogeneous medium 
typically follows well defmed laws in relation to the optical characteristics of each of the constituents taken 

3D separately. These laws are expressed in the equations for Beer-Lambert's law. The light 135 which passes through 
the finger 130 is incident upon the optical detector 140. The optical detector 140 generates an electrical signal 
proportbnal to the overall intensity of the light 135. 

Although the light 135 typically has different intensities at different wavelengths, the optical detector 140 
generates an electrical signal which is proportionate to the area contained under the spectral response curve of the 

35 light 135 within the optical band detected by the detector 140. That is, the optical detector 140 recehres light 
havmg different intensities et different wavelengths. The detected wavelengths are restricted over a band of 
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approximately 850 nm to 1.700 nm duo to th. characteri«i« of the daUctor 140. that if mtenaity » p|,ttt.d 
as a f uncton of wavelength to obtain a spectral response curve, the area under the .peetr.1 naponse curve wM be 
md.catlve of the averege optical radiation intensity incident upon the detector 140. Th«^ the electrical signal 
produced by the detector 140 is proponional to the overall (i.e.. average) intensity of the Kght 135. 

The etectrieai signal output by the optical detector 140 is amplified through the pre-amplifier 150 and is 
than passad for filtering to the bandpass filter 160. The bandpass filter 160 serves to eliminate high end low 
frequency noise signals which are extraneous to determining the blood glucose level within the patient An enatog 
agnal output from th, bandpass filter 160 is converted to . digitel signel within the analog-to-digitel converter 170 
ThB d.grtel signal is then transmitted to the digital signal processor 180 for procesting. The method employed by 
the digitel signal processor 180 to determine blood glucose level from the digitel signal provided by the analerto- 
digrtal converter 170 wiB be described in greater deteO below wi«, reference to Rg„^, 3-9. FtaaBy. a value 
mdicatree of the blood ghicose level is output from the digital signal processor 180. ' 

It should be noted here that the signal conditioning apparatus and method used to reduce noise in the sighal 
provdad to the digital signel processor (DSP) 180 advantageous^ provides a very high signaMo-noise ratio (on the 
order of 90-1 00 dB. in one preferred embodiment). Such e processing system is described in International Publication 
No. WO 96/12435 entitled SIGNAL PROCESSING APPARATUS. Particutar attention is drawn to Figures 11 11A 
and 1Z and the corresponding description of the Intermitianal Publication, which detail the structure and method used 
to obte« a signal-tOHioise retio sufficient for the application of the present invention. It dmuU further be understood 
that ahhough the International Pubncation is directed to an apparatus and method for obtaining blood oxygenation 
levels, the same signal processing steps are performed in accordance with the present invention until the 20 KHz 
data signal is output to the DSP as depicted in Figure 12 of the International Pubfication. ,. 

Figure 2 is e data flow diagram which very generaUy represents the overall method of signal enelysis 
amployed in accordance with the present invention to obtain the blood glucose level of a patient through noninvasive 
.pt.c.1 momtoring. First as indicated within an activity block 200, a patient's finger 130 is Oiumlnated by light 
wtach B filtered in a known pattern by the filter 120. The optical signal which passei through the finger 130 is 
converted to an electrical signal by the detector 140. es indicated within the activity blo^ic 203. 

In a first analysis step, electrical noise due to the electrical components of the system 100 and ambient 
Bght incident upon the detector 140 are removed es represented by. an activity block ilO. Thus, for example. 60 
Hz electrxal no.se end room Kght are removed in this step. Since ambient light and electrical noise is lineerly mbced 
w«h the desired signal, this noise can be removed by e finear subtraction at this stage in th. signal processing 
procedure. 

One. electrical noise is removed from the etectrieai signal indicative of the light 135 pessed through the 
filter 120 end the finger 130. a second anahr«s step Is perfom,ed. as represented by an activity block 220 whereby 
the electrical signal is processed to demodulate the attem.ation cheracteristics due to the fyter 120. It should be 
understood that the filter charactaristics are modulated together with the fmger characteh«ic so that removal of 
.he charecteristicsdue to the filter 120 is advantageously accompBshed by meens of demodutotion rather than Bnaar 
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subtraction. In this way the portion of the electrical signal due simply to the attenuation characteristics of the finger 
130 is isolated. 

Once the electrical signal has been processed to isolate the information solely relating to the optical 
attenuation characteristics of the patient's finger 130, those portions of the electrical signal which are due to the 
5 attenuathre characteristics of the finger constituents other than glucose and water are extracted in a third analysis 
step as indicated within an activity block 230. In the preferred embodiment, a first sub-step of this third analysis 
step invohres expressing fluids from the fingertip to teave substantnliy non-fluid elements of the f ingertqi whereby 
those finger constituents other than blood are extracted from the electrical s^naL Another sub-step in the third 
aneiysis step invoh/es extracting from the remaining signal representing blood those portions of the electrical signal 

10 which are due to the attenuative characteristics of blood constituents other than glucose and water. Both sub-steps 
are advantageously accomplished by means of demodulation to remove the signal characteristics which are not due 
to the glucose or water within the blood. 

Finally, as represented by an acthrity block 240, the ratio of the concentration of glucose to the 
concentration of water within the patient's finger 130 is taken to extract signal characteristics due to path length 

15 through the finger 130. In this manner an absolute glucose concentration is obtained. This concentration level can 
be output as a value indicatwe of the blood glucose level of the patient. 

Figure 3 is a data flow diagram which details the method used during run-time to obtain the blood ghicose 
level of a patient based upon detected optical signals. It should be understood, however, that prior to run-time, 
initialization and calibration routines other than the housekeeping and self testing procedures of block 305 are 

20 performed. 

PRE-RUN TIME INITIALIZATION 

The mitiaraation and calibration are performed at the factory or other time prior to use for patient monitor. 
In general, a blood constituent matrix b constructed and a filter characteristics matrix is also constructed, as 
described in greater detail below with reference to Figure 7. The blood constituent matrix represents the absorption 

25 characteristics of each of the main constituents within the blood over various wavelengths of light. The filter 
characteristics matrix represents the absorption characteristics of the dichroic filter 120 at different portions of the 
filter 120 and for various wavelengths of light. The filter characteristics matrix is used in order to extract portions 
of the electrical signal generated by the detector 140 which are due simply to the optical attenuation caused by the 
fQter 120. The blood constituent matrix is used for cross-correlation of light passed through tissue in order to 

30 calculate the concentration of the different constituents including glucose, water, hemdjBlobin, etc. 

The bbiod constituent matrix and the filter characteristic matrix are both two-diinensional matrices. The 
blood constituent matrix inchides one cohimn for each constituent of blood considered and one row for each 
wavelength of Gght which is measured. In a preferred embodiment, the blood constituent matrix inchides five 
columns and sixteen rows when five blood constituents are analyzed and sixteen wavelengths are used. The filter 

35 characteristic matrix includes one column for each wavelength of light which is measured and one row for each 
rotational position of the filter 120, at which an anaiog-to-digital sampling is performed. Thus, in one embodiment. 
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the filter characteristic, matrix includes 16 column, and 256 row, when 16 wavelength, are u»d and 256 rotational 
posmon, of the filter 120 are defined. It should be understood here that it i, not neces,ary that 18 different 
wavelength, be used; however, the use of additional wavelengths is particulart, advantageous for increa-ng the 
signaMo.noi»i ratio. Since about half of the incident light is transmitted through the filter at each position of the 
5 filter, the same wavelength b detected multiple times (although in a unique combination with other wavelengths each 
tene) » that the overall «gnal intenrity i, from 10 to 100 time, the intensity of any single wavetength and much 
h«her th«, the m»,e floor. This I. commonly referred to as Felgate's advantage. In this manner the spectral 
response of the entire filter 120 over the expected meamed wavelengths i, completely characteriied. The method 
employed to construct the blood constituent matrix and the fflter characteristic, matrix is described in dotal below 
10 with reference to Figure 7. . .'• 

tiUW-TIMg PBnnFsgin iR 

Once initialization end calibration have been performed, the system is ready for run time use. A, depicted 
.n FQure 3. the start of run-time processing is represented in a begin block 300. First! housekeepinB and self -testing 
procedures are performed, a, represented in an acthrfty block 305. Briefly, housekeeping and self testing involves 
hoot operation, ami ensuring proper operation of the btood glucose moratoring system 100, and to provide more 
accurate monitoriog of the blood ghicose concmrtration. For example, the instrument first determines if there is a 
auffic-nt sigmri intensity to take an accurate reading, is there a heartbeet. is the petient following instructions by 
pressmg and releasing the finger 130. etc. After housekeeping end self testing is completed, the light source 110 
« activated to transmit light 115 through the filter 120, as represented in en acthrity btock 310. InitiaBy. the light 
source 1 10 is actuated while the patient's finger 1 30 is not interposed between the filter' 120 and the detector 140 
Thus, the light which is detected by the detector 140 represents a baseline light intensity (U which can be used 
as a test to in«.re that a bulb which is too dim or too bright is not inserted as a replacement bulb for example. 

Once the faiitial baseline light intensity constant has been determined, the patient inserts the finger 130 
as mdmated in an activity block 312. so that measurements of the blood glucose level within the patient's finger 
130 may be tricen. A, described above, when the patient', finger 130 is inserted between the filter 120 and the 
detector 140. Sght 115 from the source 110 passe, through the filter 120 end the finger 130 to be detected a, 
Bght 135 incident upon the detector 140. 

As indicated within an activity block 315, the light which is incident upon Ihe detector 140 is comrerted 
to en electrical signal and this signal is ampBfied in the pre-amp 150, filtered with the band pass filter 160 md 
s»npled by the enalog-to-digitat converter 170. Since the filter 120 is rotating (at approximetely 78.125 revohitions 
per second in one actual embodiment, although other rotational rates could be advantageous as called for by the 
partmular appication). samples of the electrical sigm,l output by the detector 140 ere indicative of the light intensity 
otected at v«iou. rotational positions of the filter 120. In one advantageous embodki»„t. one complete rotation 
(..a, 360 ) of the filter 120 correspomi. to 512 digitel samples. That is. 512 ,ampki, are taken within the period 
correspondingto one revotation of the filter 120. Thus, f or exmnpto. if the f Uter 120™tetes et 78.125 revolutions 
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per second, then 512 samples will be taken within approximately 1/78th of a second, so that the sampling rate of 
the analog-to-digital converter 170 wlQ be approximately 40,000 samples per second, i 

As further described below, the filter 120 constructed in accordance with the present invention inchides 
redundant regions within an entire revolution. Spectficafly, the filter 120 is symmetricaily layered so that the first 
5 half-revohition of the filter provides a mirror of the signal of the second half-revolution of the filter 120. That is 
to say, as depicted in Figure 10, the fitter is formed in a wedge shape so that the thickness in one direction is 
constant and the thickness m the perpendicular direction increases fmearly. Thus, the second half-revolution of the 
fitter 120 is redundant. For this reason, digital samples taken for one-half of the revolution of the fOter 120 could 
be discarded so that in each rotation of the filter 120 there are 256 samples used for purposes of digital signal 

10 processing rather than 512 samples m the embodiment described above. Alternatively, all 512 samples can be used 
for processing by averaging corresponding values. In yet an alternative embodiment, the redundant half of the f Bter 
may be used for filter and source calibration. Each of the 256 samples (if only half are used) represents a different 
portion of the fitter 120 having different optica! transmission characteristics. 

Once the signal output by the detector 140 has been ssmpled, linearly added determ'mistic and stochastic 

IS electrical noise inherent withm the blood ghicose monitoring system is extracted (le., linearly, subtracted), as indicated 
by an acthrity block 320. The method of subtracting this noise depends upon whether the noise is deterministic or 
stochastic. If the noise is deterministic, then this can be modeled with the appropriate phase and subtracted. If 
the noise is stochastic, then the noise can be averaged towards zero. As described above,!the filter 120 is specially 
designed to include an opaque strip (i.e., the brass strip 122). The digital signal processor 180 detects when the 

20 opaque strip 122 of the filter 120 is interposed between the light 115 and the detector 140 by monitoring the 
intensity output from the detector. This intensity will be effecthrely zero when the light is blocked by the opaque 
strip 122. Since the opaque strip 122 blocks substantially all of the optical radiation transmitted from the source 
110, any signal output from the optical detector 140 when the Bght is blocked (e.g., from ambient light, thermal 
effects, etc.), will be interpreted as electrical noise which is not due to either the spectral absorption characteristics 

25 of the finger 130 or the spectral absorption characteristics of the filter 120. Thus, the digital signal processor 180 
interprets the signal present at the output of the optical detector 140 when the brass strip 122 is interposed 
between the light source 1 10 and the optical detector 140 as stochastic noise which is subsequently subtracted from 
all signals output from the optical detector 140. In one embodiment, this is simply accomplished by subtracting the 
digital value corresponding to the detected noise level from each of the digital values corresponding to the detected 

30 signal samples obtained within the acthrity block 315. Alternatively, a shutter mechanism could be interposed within 

the Sght path, or the lamp 110 could be turned off momentarily to provide the same effect. In th» manner, the 

> 

electrical noise inherent within the blood glucose monitoring system 100 is removed so that those electrical signals 
due to the optical absorption characteristics of the filter 120 and of the finger 130 are considered in the further 
processing steps. 

35 Once the stochastic noise inherent within the blood glucose monitoring system 100 has been extracted by 

averaging to zero and the deterministic noise is subtracted by phase modeling the noise or averaging the deterministic 
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noise to zero, control of the method passes from the activity block 320 to an activity btock 321 Within the activity 
block 323 the signal is divided by 1. to nomialize the signaL The nom.eittd signel is subsequently processed within 
an act«ity block 325 to construct a signal intensity matrix, or vector, from the sampleivalues obtained witiun the 
actwrty block 315 (taking into consideration the subtraction of the electrical noise performed within the activity block 
320 end the signal normalization performed in the acthrity block 323). The signal iminsity matrix is a one column 
matnx (sometimes referred to as a vector} incbding 256 signal intensity welues (i.e.. one value for each sampled 
rotational position of ti» filter 120). Thus, the signal intensity vector Is obtiiined by direct measurement of tt,e 
optical signal which passes tiirough boti. the filter 120 and the finger 130 end is deti»ted by the optical detector 
140. Of course, the veh.es used to fom, ti,e sigmil intensity vector are taken from tii. .npitude of the tignah 
output from the detector 140 after subtraction of the «,isa from each sample. DasiBnatin, each rotational position • 
of a fdter 120 which is sampled by the analog-to-digital converter 170 by the symbol «. ti,an 0, wl correspond 
to the fust rotational position of the filter 120, <p, wffl correspond to the second rotational position of tiie fifter 120 
to 0«, which corresponds to the last rotational poshion of the fBter 120 before 0, is taken again. Using thii 
notatmn. I„ corresponds to the intensity of light detected by the optical detector 140 when the filter 120 is in the 
first rotational position 0„ I,, corresponds to tiie intensity of light detected by the detector 140 when the fiher 
120 is in tiie second rotational position 0, etc, Thus, the signal intensity matrix comprises a single column matrix 
havmg 256 digital values from I,, to W which correspond to ti,e opticel intensities detected at each of the 
rotational positions of ti» filter 120. In om, embodiment, ti.e intensity vah.es for several revolution, are averaged 
to form the signal intensity matrix. 

Subsequent to the construction of the signal intensity vector withm the actwityfhtock 325. error checking, 
or testing, is performed within an activity block 327. This test is periodically performed to insure that a void finger' 
sampte is being monitored. Thus, at each stage of transformation and fihering. the output is examined for en 
expected renga of values to corresponding to the expected renge of properties for a human finger. If the output 
values are found to be off-scate. then a system error occurs to indicate that the monitored sample is not valid 

Once error testing is performed on the signal intensity vector, hereinafter designated as 1(0). and the fiher 
characteristics matrix. heraim.fter designated as F(0.^). has been obtained prior to run-time. as indicated above ami 
represented as a data input in a block 333. tiie signal intensity matiix together with the filter char«:teristfcs matrix 
may be used to obtain a matrix indicative only of the optical absorpttan charactaristh:, of the finger 130 as 
represented in actWity blocks 330, 331. That is, since the overall optical absorption characteristics of both the f'iher 
120 and the finger 130 are known as measured within the signal intensity matrix. I(^). and the optical absorption 
Characteristics of the filter 120 are known as represented by the filter characteristics matrix, ^<p^) the optical 
abaorption of the light 115 due to the characteristii» of the finger 130 may be determined by removing the optical 
absorption due to the fdter from the over.il absorption of ti,e filter 120 and the finger 1^0. This b accomplished 
by fust takhtg ti» mvarse transform of ti.e filter matrh, as represented in tite activi,y>ock 331. end suhsequentiy 
n«.ltq,lying the eignal inti,nsity vector by ti,e inverse filter matrix, as represented in the activity block 330 
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Figure 5 graphically depicts the relationship between the optical absorption due to the finger 130, the 
optical absorption due to the filter 120, and the total optical absorption due to both the finger 130 and the filter 
120. As depicted in Rgure 5, transmission of 6ght over wavelengths from 850 nanometers to 1,350 nanometers 
through the finger 130 results m a spectrum function designated as TM) wherein the transmission of light through 
the finger 130 is plotted as a function of the wavelength. In similar fashion, the transmission of light through a 
selected rotational position (e.g., when ^ - 0, corresponding to 0*") of a filter 120 is plotted as a function of 
wavelength and is des'^nated by the function F(0^) in Figure 5. Fmally, the combination, or convohition, of the 
optical absorption due to the finger 130 and the fitter 120 is shown over the same wavelengths and represented 
In Figure 5 by the function \l<p). To obtain from the finger transmission function TU) and the filter transmbston 
function the optical transmission percentage at any wavelength along the functions TM) and F(0^) is 

multiplied to obtain \{<pl Thus, for example, at a wavelength of 1,050, the transmission of light through the finger 
130 is approximately 0.24%, while the optical transmission through the fiher at the same wavelength for <p ^ Q'' 
is approximately 80% so that the total optical transmission through both the finger l3b and the filter 120 will be 
approximately 0.24% times 80%, or a total of 0.192%, as indicated at the wavelength 1,iD50 and the function 1(0). 

The functions 1(0) and f{0M may be represented by the signal intensity and filter characteristic mauices, 
respecthrely. Thus, since 

/(<!)) = F((j),A) X T(X) (3) 

and represents a one-column matrix (vector) containing an intensity value for each rotational position value 0, 
while represents a two dimensional matrix containing a filter transmission coefficient vahie for each value 
of 0 and each value of A (see Hgure 4DI, then the function JiA). representative of optical transmission through the- 
finger 130, may be represented as a one column matrix having values for each of the various wavelength values, 
A. K ' 

In accordance vvith one embodiment of the present invention, 16 wavelengths are selected over the range 
of 850 nanometers to 1,400 nanometers for purposes of characterizing the spectral characteristics of the finger 130 
as well as the filter 120. Specifically, in the preferred embodiment, the monitored wavelengths are 850, 880, 910, 
940, 970, 1000, 1030, 1060, 1090, 1120, 1150, 1200, 1250, 1300, 1350, and 1400 nanometers. 

The matrix form of equation (3) above is shown below: 

As shown in Equation (4), the signal intensity matrix 1(0) is equal to the product of the two dimensional filter 
characteristic matrix, H^Al and the single column finger characteristic matrix HA), In this equation, two of the 
mauices are given (i.e., 1(0) and HpMl Thus, the third matrix, llA), which represents the optical transmission 
characteristics of the finger 130 for the 16 selected wavelengths between 850 nanometers and 1,400 nanometers, 
may be obtained by simply multiplying the inverse of the filter characteristic matrix, designated as F^0M), by the 
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signal intensity matrix. using conventional matrix inversion and multiplication technigues, as shown below. 
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Thus, as indicated in an activity block 331. the inverse transform is taken of the filter characteristic matrix. FHipM 
and then this inverse matrix is multiplied by the signal intensity matrix. 1(0). within the acth,rty block 330 to obtaiii 
the frequency response of the finger 130 as expressed by the finger characteristic matrix, or transmission vector 

Once the transmission vector (or. as aftamately referred to herein, the optical frequency response matrix) 
TM). has been obtained as indicated within the acthrity block 330. the digital vahias stored for each of the 16 
selected wavelengths within the transmission vector are converted from a substantiallV logarithmic function to a 
Imear function, as represented in a subroutine block 335. The method for comrerting the substantially logarithmic 
values stored within the transmission vector. TM). to linear vah.es wiU be described in greater detail below with 
reference to Figure 8. 

lINEAHIZflTinw OF THP TBflNSM|SSlQN VKTOR : 

The uansmission vector. TM). Is linearized to prevent cross modulation of interfering signals. By linearizing 
the transmission vector, characteristics of the detected optical signal which are not due to the optical properties of 
tl» glucose Within the bloodstream can be lineariy subtracted. Thus, the sequence of analysis and Imeerization u«.d 
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in accordance with the present invention is very important to prevent cross modulation of interfering signals with 
the desired glucose signal 

The ultimate purpose of Gnearlzlng the transmission vector, HA), is to e&ninate path length dependent 
errors. As discussed above, the transmissran intensity of fight through a medium (e.g«, such as the finger 130) has 
5 a substantially exponential relationship with the path length of the light transmitted through the medium. Thus, the 
different transmission characteristics vary non-linearly with the corresponding path lengths so that taking an average 
transmission intensity value over each of the possible transmission path lengths will not result in a transmission 
intensity value which is independent of path length. This is because, due to the non linear character of the relation 
between the path lengths and the transmission intensity values, some path lengths are weighted more than others 

10 so that errors along the heavily weighted paths become more pronounced in the final blood glucose concentration 
output value. By iineariznig the optical frequency response matrix, the path lengths become linearly related for each 
wavelength so that the path lengths can be subtracted from the set of linear equations. Once the values within the 
transmission vector, TM), have been finearized, this linearized matrix is designated by the new name DM) herein, 
which represents the optical density. 

15 Once the values within the transmissmn vector have been linearized as indicated within the subroutine block 

335 (Figure 3), control of the method passes to an activity block 340, wherein, in accordance with one aspect of 
the present invention, significant pressure (i.e., somewhat above the patient's blood pressure) is applied to the finger 
tip 130 to express fluids from the finger tip 130. It should be understood that this step in the monitoring process 
is part of a long loop method which is not performed along with the other signal processing steps. This is because 

20 static conditions are required each time the simultaneous equations derived from the optical density vector are solved. 
Thus« pressure is applied to the finger for a long enough duration to derive a solution to one or more simultaneous 
equations, and then the finger is released for another period of time sufficient to obtain accurate solutions for 
additional simultaneous equations. Thus, since the application of pressure to the finger tip 130 is performed at a 
slower rate than the signal processing steps, the finger tip is, for practical purposes, in a static condition during each 

25 iteration of the signal processing method. Fluid is expressed from the finger tip 130 to obtain measurements of the 
optical characteristics of the finger not due to blood. This allows removal of portions of the electrical signal 
generated by the detector 140 corresponding to optical characteristics of the f'mger not due to the blood within the 
finger tip 130. Briefly, the expression of fluid within the finger tip 130 involves appiying physiological pressure or 
other means to the patient's finger 130 and releasing the pressure so as to cause substantiaDy all of the blood to 

30 evacuate the finger 130. As is known hi the art, the finger 130 contains not only blood, but also bone, flesh, skin 
and nail, all of which absorb optical radiation, and therefore, contribute to the absorption characteristics defined 
within the transmission vector. However, for purposes of the present invention, it is not useful to mclude the optical 
absorption characteristics of bone, flesh or other non-blood constituents of the finger 130 for calculating blood 
glucose; these finger constituents may be considered as artifacts which should be extracted from the data used to 

35 determine blood glucose levels. 
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In order to remove these artifacts due to non-blood constituents of the finger 130, significant pressara is 
applied to the finger tip 130. When the finger 130 is physiologically altered so as to cause blood to flow in and 
out of the finger tip 130. the optical t»n«iussion characteristics of the finger tip ISO'when the blood is evacuated 
from the fmger tip 130 diNers from the optical transmission characteristics of the fingertip 130 when blood flows 
into the fingertip 130. The difference is effecthrely the optical transmission characte^cs of the blood, and other 
fhiid. which was evacuated. Stated another way. the optical characteristics of the fingartip 130 which is fufl of 
blood is a combination of the optical transmission characteristics of the fingertip 130 when the blood is evacuated 
and the optical transmission characteristics of the evacuated blood. Thus, by subtracting the signal produced by the 
optical detector 140 when blood is evacuated from the finger from the signal produced by the optical detector 140 
when the finger is full of blood, the opticel transmission characteristics of the evacuated blood may be obtaimid 
It will be umlerstood by this description that the blood need not be completely evacuated. In this manner, artifacts 
due to bone, flesh, nay and skin ef the fingertip 130 ere removed, or extracted, from the digitally processed signal 
ge«,rated by the opth:al detector 140. It should further be noted that other methods to express blood from the 
finger tip 130 may inchide, for example, raising the patienfs arm in the eir so that gravity causes a decrease in the 
IS flow of blood in the finger 130. 

Figure 12 depicts an exemplary embodiment of a device which could be usedtd^provide for expression of 
fyds from the fingertip 130 during spectral analysis in accordance with the present invention. As shown in Hgore 
12. a hand-held glucose monitor 1300 includes the light source 110 , the lens 117. the dichroic filter 120, and the 
detector 140 (aU shown in phantom). The fingertip 130 is placed between a caliper pressure arm 1310 which may 
be used to epply pressure to the fingertip 130. A display 1320, which displays glucose concentration, is also shm»n 
in Figure 12. 

In a particularly adventegeous embodiment, a pronounced perturbetion of the fingertip 130 may be 
performed in order to obtain additional information regarding the optical characteristic, of the fingertip 130 and the " 
blood therein. In such an embodinent. the fingertip could be placed within en inflatable jacket which is mflated to 
e relatively high pressure to squeeze most of the blood out of the fingertip 130. Thus, nearly 100% modutetion of 
the signal due to the blood flow in the fingertip 130 is achieved in this manner. This high modototion increase, the 
.ignal-to-noise ratio of the electrical signal due to the blood in the fingertip 130. Of course, it should be noted thet 
such hQh modulation may change the path length through the fingertip medium (since the thickmiss of the finger 
mey be significantly ahered). Such a change in path length may require compensatioi,.; 

For this reason, b, one embodiment of the bivention. a second water coef fici^ vakie is included in the 
blood constituent, matrbi. hereinafter designated as AM.), where r designetes the nuhier ef blood constituents 
which ere used to form the btood constituent matrix. This second water coefficient value may be used to determme 
the change in the "depth- of water through the fingertip when there is significant modulation of the fingertip Thus 
a f.st water depth is calculated using e first coefficient at one time when the finger is not «,u.ezed. white a second 
weur depth is calculated using the second coefficient at another time when the finger is squeezed. The change of 
the depth of the water m the fingertip is substantially the same as the change in path length through the finger 
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so that the change in water depth may be taken to be representative of the change in path length. Thus, by using 
first and second coefficients for water in the finger 130 at selected intervals, a ratio may be obtained which is 
indicathre of the ratio between the path length through the finger 130 with modulation and the path length through 
the finger 130 without modulation. 

Once the signal features due to the non-blood artifacts within the fingertip 130 are extracted, those optical 
characteristics which are due to the biood constituents other than glucose and water are extracted from the optical 
signal, as represented in an activity block 345. The blood glucose concentration level within the water of the blood 
stream is then determined, as further explained. 

As indicated within the activity block 345, a concentration matrix is derived from the linearized finger 
frequency response characteristic matrix, D{M (after extraction of nonblood signal features) and the blood 
constituent matrix, AU.r). The blood constituent matruc is provided as a data input to the 'brocess performed within 
the activity block 345, as indicated within the bhick 347. As is well known in the art, and discussed briefly above, 
when several substances are combined to form a medium with optical absorption, the overall optical absorption and 
transmission characteristics of the combined medium are typically describable in terms of the optical characteristics 
of each of the separate constituents of which the medium is composed. Thus, for example, if a medium is composed 
of water, oil, and alcohol (which are only used here for illustrative purposes, since, in practice, these elements are 
insoluble), the optical transmission through the entire medium may be described, in accordance with Beer-Lambert's 
law, as proportional to the sum of the exponents of the optical absorption characteristics of the first constituent 
fte., water) times a concentration of the first constituent times the path length through the medium; plus the optical 
absorption coefficient of the second constituent (Le., oil) times the concentration of the second constituent times 
the path lengths through the medhim plus the absorption coefficient of the third constituent n.e., alcohol) times the 
concentration of the third constituent times the path length through the medium. That is,. 



Since the optical absorption coefficient for each material typically varies as a function of wave length, a 
series of equations may be derhred to describe the overall optical absorption characteristics of the entire medium as 
a function of the constKuents of that medium. Thus, 

Dividing both sides by Iq and taking the logarithm of each side provides the following equation: 
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for , ha f«, w.v.lB„Bth>»,. for exampte. If the optical trans™«io^ 

Boer-Lamberfs law. than a first loaarithn. would result m a set of linear equations as shown above. From these 
equ-uon,. aimple matrix algebra could be u«.d to obtain the concentration of olucose h. blood. However as will 
be dBcussed in flraater detafl below, the optical characteristics of multiple constituent, within the bloodsueam do 
not result in a set of linear equations after takino , sinBia logarithm. Thus, the present Invemion contemplate, a 
erent approach which includes the steps of taldno a near-log of the values within the transmission v«:tor. 
foUowed by two true logarithms, or another near-log in an alternat«e embedment, to obtein more precise linearisation 
of the transmission vector. It should be understood that preWnary linearization and other correction sttps are used 
n. co.*.nctio„ with the above steps to insure that the conditions of Beer-Lambert's law are satisfied before either 
logar«hm Is taken. Furthermore, it may not be necessary to perform the second logarithm if a sufficiently linear 
outcome is observed after the \nx logarithm is taken. 

The double logarithm process is used in eccordance with the teachings of the present invention to arrive 
Ot a senes of linear equations which may be described in matrix fom. This Inearization takes place wrthin the 
subroutme block 335 and is described in greater detaS below with reference to figure 8 

The transmission vector. TM). is now designated as a smgte cotamn taearized finger spectrum matrix. OM) 
Thus, for each of the ,6 selected wavelengths of light, there exists a digital value corresponding to the linearized 
opt«al frequency response matrix, m. That is. DM) comprises a first value D. indicating the linearized absorption 
.o«o.nt for the finger f, ,30 at wavelength ... D« indicating the Bnearized attenuetion coefficient for the finger 
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AUrfh and the concentration times path length matrix, hereinafter designated as C(r)PL This representation is shown 
in matrix form as: 
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where PL represents the path length and c,, Cj. c, represent the concentration of the first blood constituent, the 
second blood constituent, and the last biood constituent (designated by "r"). In accordance with well known 
matrix algebra techniques, the matrix representing the concentration times path length matrix, C(r)PU may be 
determmed by taking the product of the linearized optical frequency response matrix, DU), and the inverse of the 
blood constituent matru, AMrri, as depicted below. 
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As indicated within the activity block 345, the inverse transform, A^Urr), of the blood constituent matrix 
together with the linearized optical frequency response matrix may be used to obtain a matrix indicative of the 
concentration of the different constituents within the patient's blood. 

Once a matrfac indicathre of the concentration of the different blood constituents has been derived as 
mdicated within the acthrity block 345, the ratio of glucose concentration thnes path length through the finger 130 
to water concentration times path length through the finger 130 is calculated as indicated within an activity block 
350. This ratio is the concentration of glucose in water which is the same as the glucose concentration within the 
blood stream. The above described method is performed over several iterations so that several glucose concentration 
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values are obtained. To obtain an average glucose concentration value, least squares analysis, which I, well known 
to those skied in the art is used to ptot a line indicative of blood glucose concentration as indicated within an 
actnnty block 355. This value is output to a dbplay (see Figure 121 which may be read by the patient or by the 
operator of the blood ghicose n»nitor system 100. Control of the method then passes from the acth,ity block 355 
5 to an end bbick 360 wherein the method of the present invention is completed. 

PPRIVATIQN OF THF FHTFR rHftRACTCRlsnn Awn rloqd cauvmif ftfj fff^r^^^«> 

Figures 4A4D, together with Figure 7, iflustrate in greater detail the method employed in accordonce with 
the present invention to perform pre run-time initialization and calibration of the blood ghicose monitoring system 100 
The rmtmlization and calibration routine is illustrated in Figure 7 ami starts with a begin block 800. " Control passes 

10 from the begin block 800 to an activity block 810. 

Activity blocks 810.825 Olustrate the method employed in accordance with the present invention to obtain 
the blood constituent matrix. AM.r). white activity blocks 830.845 illustrate the method used to obtain the filter 
characteristics matrix. It should be understood, however, that although the method depicted in Figure 7 is 
represented as ti,ough the blood constituent matrix is constructed folowed by the filter characteristics matrix, in 

15 Pfactice these represent independent procedures which may be run in paralleL 

As discussed above, different substances typ»ally have different spectral charecferistics. That is. different 
substances absorb light more or tess across different wavelengths of light Thus, when the optical transmission of 
light .s plotted versus wavelengths for a given substance, a pattern, which is sometimes refwred to as ti,e spectral 
SHinature of a substance, is observed. This signature defines the spectral characteristics of the substance over the 

20 various wavelengths of light. 

To form the btood constituents matrix, AM,rl, the main constituents of blood are separated into individual 
substances, and the spectral characteristics of these substances are plotted versus wavelength. Separation of the 
mdivaiuel substances may be performed by an actual physical separation of the constituents in vitro or by means 
of an « clinical test wherein the concentrations of the individual substances are carefully controHed to provide 
a reference for monitoring the concentrations of the various blood consthuents. In the,pre«,nt embodiment since 
16 wavelengths between 850 nanometers and 1,400 nanometers are used to characterizttim optical characteristics 
of the pat«nt's finger 130. the same 16 wavelengths are used in order to construct the btood constituents matrix. 
Thus, for each substance or main constrtuent within the blood, a measurement is taken to obtain the absorption at 
each of the 16 wavelengths. 

Each of these 16 velues. which are referred to as spectral absorption coefficients, is used to define a row 
of the blood constitiients matrix, AM.r|. Likewise, each column of the bteod constituenu matrb. will correspond to 
a particular blood constituent so that there is the same number of colomns in the blood constitiients matrix as tiiere 
are blood constituents. Thus, if there are V constituems of btood and V wavetengtiis over which the optical 
character,st.cs of these constituents are to be determined. ti» bteod constituents matrix comprise, an "n" row by 
35 "r" column matrix. 
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As indicated within the activity block 810, for the cases of biood constituents that can be physically 
separated, the first constituent (e.g., water) is illuminated as a separate substance over the n wavelengths of fight 
(i^^ in this embodiment, 16 wavelengths) to obtain the spectral absorption coefficients of water over these end 
wavelengths. For example, a calibration sample could be used in this case. Then, as indicated in an activity block 
5 815, the second constituent of blood is ilhiminated with the n wavelengths of light to obtain the spectral absorption 
coefficients of the second constituent (e.g., oxyhemoglobm). This is repeated for each of the main blood constituents 
until the "rth" blood constituent (eg., glucose in water) is Illuminated over the n wevelengths to obtain the spectral 
absorption coefficients of the last main blood constituent 

In a preferred embodiment, certain constituents, such as water and glucose dissoNed in water, are measured 

10 in vitro, while other blood constituents, such as oxyhemoglobin and deoxyhemoglobin, are preferably measured in vivo 
due to the difficulty of obtaining accurate in vitro measurements for these constituents. Furthermore, it should be 
noted that m a non-imaging system, such as typically used to perform in vivo measurements, the source and the 
detector should be diffuse. The in vhro and in vitro measurements are used as an approximation, and an iterathre 
process is used with the approximate values until the Imearization process described above converges. In this 

15 manner, n (e.g., 16 in the present embodiment) absorption coefficient values are obtained for each of the r 
constituents. In accordance with the present invention, one preferred embodiment the spectral characteristics for 
ghiCGse, water, oxyhemoglobin and deoxyhemoglobin as separate substances are determined. In addition, a fifth 
"constituent" row within the blood constituent matrix is defined by absorption coefficients of blood due to scattering. 
A further sixth constituent row may be added where the sixth constituent is water again. This addittonai constituent 

20 may be used to double check the peth lengths obtained at different times during, for example, application of 
significant physiological pressure to the fleshy medium. 

As is well known in the art, the transmission of optical radiation through a medium typically mvohres the 
scattering and reflection of light waves through the medium. This scattering effecthrely increases the path length 
through the medium. Since some of the effects due to scattering are statistically well defined, the extinction effects 

25 due to scattering may be treated as an additional constituent within the bloodstream. The ''absorption" coefficients 
for the scattering "constituent" may be determined empirically by transmitting light through an optical medium of 
which the concentrations of each of the constituents is already known. 

Once ell the main blood constituents have been defined over the wavelengths of interest for eech of the 
16 wavelengths from 850 nanometers to 1,400 nanometers, the blood constituents matrix is formed from the 

30 spectral absorption coefficients obtained for each blood constituent. Thus, the first cohmm of the blood constituent 
maUix has each of the ebsorption coefficients for wavelengths /t1, M to ^n, for the first blood constituent (e.g., 
water). The second column of the blood constituent matrix has the absorption coeffuiieht vahies for each of the 
wavelengths yl1, At to yln, for the second blood constituent (e.g., oxyhemoglobin). Aiid so on until the last column 
of the blood constituent matrix comprises the spectral absorption coefficients for the last main blood constituent at 

35 the selected wavelengths Ah At ^, An. 
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The activity blocks 830-845, together with Figures 4A4D, illustrate the method used in accordance with 
the present invention to construct the filter characteristics matrix. As discussed above, the filter 120 reflecta and 
transmits optical radiation in different proportions for different wavelengths at different places on the filter disk 120. 
This b dearly illustrated in Figure 4A4C. wherein Figure 4A represents the optical transmission of Gght at a 
wmrelength of 850 nanometers ptotted versus each of a possible 256 disk rotational poshions. A, shown in Figure 
4A, when the disk 120 is in the initial starting position fte.. p - O). the transmission of fight at 850 nanometers 
is approximately 10* through the filter 120. whfle when the disk 120 is rotated so that 0-32, the optical 
transmission of light at 850 nanometers through the filter 120 is approximately 25%. Again, between the disk 
rotational positions of 0 - 128 to « - 160. the transmisston of Gght at 850 nanoriwer, wavelength throuoh the 
fflter 120 is approximately 75%. Thus, the optical transmission for A - 850 nanometer, is entirely characteriied 
over 256 rotational positions of the disk fiher 120. as depicted in Figure 4A. 

Figure 48 depicts the optical transmission characteristics of light at 1.150 nanometers over the same 255 
rotational positions of the disk 120. Simiarly. Figure 4C depicts a plot of the optical transmission of fight at 1 350 
nanometers through the disk filter 120 at each of the 256 rotational positions of the disk 120. In one actual 
embodiment of the invention, the optical Uansmission characteristics of the filter 120 are described for 256 rotational 
positions at each of 16 wavelengths between 850 nanometers ami 1.400 nanometers. 

Thus, from these measurements, a filter characteristic metrix may be constructed, as shown in Figure 4D. 
The filter characteristic matrix designated in Rgure 40 as f{<pM mchnles 256 rows «.d 16 cohimns. Each cohimn 
of the filter characteristic matrix comprises the spectral absorption charactaristH» of the disk 120 at each of the 
20 256 rotational positions of the disk 120 for a given wavelength. 

in order to construct the fSter characteristic matrix depicted in Figure 40, the filter 120 is lUuminated at 
a first rotational position over each of the 16 wavelengths to obtain spectral absorption coefficients for each of the 
16 wavetengths, as indicated within an acthnty block 830. Once the spectrel absorption coefficients have been 
determined for the first rotational position as mdfeated within the activity bhick 830. the filter is ittuminated at a 
second rotational position (i.e. 0 - 1) over the 16 »,h,cted wavelengths to obtain spectral absorption coefficients 
for the second rotational posHion, as represented in an actnity Uock 835. This method is carried on for each of 
the possible rotational positions of the disk 120 until, es indicated within an activity block 840. the filter is 
liiummeted at the "mth." or last, rotational poshion (i.e., position 256) of the disk filter 120 over the 16 selected 
wavelengths to obtain the spectral absorption coefficients for the last rotational position. In one preferred 
«nbodiment. where a stepper motor is used, the rotational positions will be precise from revolution to revolution of 
the dnk 120. Of course, a computer disc motor with salient potes and run at a constant speed could be used 
provided that phase dithers are mmimized to less than one part in 256. 

Once spectral absorption coeffieienu have been deurmkied for eU 16 wevelengtiis of all 256 rotational 
positions of the disk 120. the fUter characteristics matrix is constructed, as indicated whhin an activity block 845 
by putting coefficients corresponding to various blood constituents in rows end wavetongti» of the transraHted Gght 
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in columns. Once the filter characteristics matrix and blood constituent matrix are constructed, the system has the 
necessary constraints for routine processing. 
THE SOLVENT EFFECT 

Although ft is typically the case that when several substances are mixed to form constituent elements of 
an optical medium, the optical transmission and absorption characteristics of the medhmi relate to the optical 
transmission and absorption characteristics of each of the constituent elements withm the medium as described by 
Beer-Lambert's law, it has been found that in certam isolated cases, this is not the case. In particular, when gbicose 
is mixed in water so as to dissolve in water, the optical characteristics of the combination of glucose and water do 
not correlate directly to the optical characteristics of glucose and water separately m accordance with the relation 
described by Beer-Lambert's law. In fact, when glucose is dissolved into water, it is found that the glucose water 
solution has a lower absorptivity (f.e., a higher transmicivtty) than water without glucose at certain wavelengths. 

In past systems which monitor blood glucose concentration by means of optical signal detection and 
processing, this solvent effect, which was heretofore unrecognized in applications mvolving blood spectroscopy, wss 
the cause of a significant amount of inaccuracy and detection of the actual blood glucose concentration. Thus, in 
accordance with the teachings of the present invention, special modifications are made to both the blood constituent 
matrix, AU,r), and the filter characteristic matrix, f{0M in order to compensate for inaiccuracies and non-linearities 
due to the soWent effect. 

Specifically, the blood constituents matrix includes ebsorption coefficients for glucose as dissolved in water, 
rather than glucose as a separate substance. The scale factor by which the absorption coefficients for glucose 
dissolved in water differs from the absorption coefficient of water is approximately equal to iog(-Iog 1^) - log(-log 
TJ, where T^ is the transmicivity of water having glucose dissolved therein, and T^ Is the transmichrity of pure 
water. 

Furthermore, the wavelengths which are to be monitored (which defines the fOter characteristics matrix) 
are selected differently than would be expected for a glucose/water mixture. This is because the solvent effect 
causes a small shifting in the frequencies of maximum absorption of the blood constituents. For example, glucose 
and water acting as separate agents would have respective wavelength absorption maxima at approximately 1070 
nanometers and 975 nanometers. However, due to the solvent effect, the water with glucose has one absorption 
maximum at approximately 960 nanometers. 

Since the absorption characteristics of pure water differ from the absorption characteristics of water having 
dissolved glucose es a function of wavelength, this difference in absorption characteristics may be used to scale 
absorption coefficient values for ghicose dissolved in water within the bloodstream. Ftgiire 11 represents the ratio 
between the optical absorption characteristics of water and water with ghicose. As shown in Figure 1 1. water 
containing dissolved glucose has substantially lower absorption than pure water at approxniately 960 nanometers, 
1145 nanometers, and 1380 nanometers, while pure water has substantially lower absorption than water containing 
dissolved glucose at approximately 1060 nanometers, and substantially the same absorptivity is observed for both 
at wavelengths of 1000 nanometers, 1100 nanometers, and 1230 nanometers. It should be noted that those values 
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Of the curve of Figure 11 within the shaded region are due to the soh,e„t effect and are ti»rafore unexpecud in 
normal spectroscopy applications. The values of the difference between the absorptivities of pure water and water 
wrth Blueo», at the peak wavelengtiis (i.e. 9B0, 1060. 1145, etc.) vary as a functton of glucose concentration 
By using the wavelengths for which the absorptivity b substantially the same for both pure water and water with 
glues, as a baaeln.. nsaGng factors for determmlng the absorption coefficients of glucose dissolved in water may 
be daterminad. Thus, the present invantion prmridas for mora accurat. measurement of the blood glucose level due 
to recognition of tha sohr«« effect, and mcraasad sensitivity if one usas 1060 nm wmlangth the minimum) 
to provide an increased difference value with tha 860 nm and 1145 nm wavelength, fte, the maximum) 
METHQP FOR iiNFftRiziWG THF npytr H FREnuFMnv i^ FS pn NffF fCmr 

Figures 6A^C together with Figure 8 depict the method used in accordance with the present invention to 
fcnearize the values within the transmission vector, TMl. As depicted in Figure 6A. the percentage transmission of 
«pt«:al radiation at a particular wavelength (each curve representing a d«eren, path length PI, P2, etc.) is plotted 
«r«.s the product of the absorption coefficient and the concentration of the optical medium through which the 
optnral radiation passes. Thus, the curve of Rgur. 6A indicates tha relation between how much Bght is transmitted 
through a g»en pathlength of a medium, and the concentration of a given substance within the medium A total 
of eight curves are sh««n in the graph of Figure 6A wherein each curve indicates a differm,, path tength {e.g. PI 
P2, ^ P8) of fight through the medium. 

As is well known in the art. some optical radiation which passes thrmigh an optical medium such as the 
fmger 130 passes immediately through the medium without substantial scattering so that tha path tongth observed 
by that portion of optical radiation is approximately the same as the thickness of the finger. In other cases however 
the optnral radiation scatters whhin the medium so that the effective path length observed by these portions of the 
optHsal rad««ion Is subnamialhr longer than the thickness of the optical medium. It has been found that the average 
path length observed through an optical medium «,ch as the finger 130 is typically three to five times the thickness 
of the finger 130. It is importent that a Inaar relation exist among the several path tengths of optical radiation 
through an optical medium such as the finger 130 in order to insure that any given path length results in tha same 
relation between the concentration of a substance within the medium and tha optical characteristics of tha medium 
« linear relation does not exist for various path lengths, a set of linear equations is not obtainable and the solution 
darned from the matrix equations wfll not be accurate. 

As .hm»n b, Figure 6A. a linear relation does not exist among the several path lengths of optical radiation 
«.r«.Bh an optical medium. Rathe,, each of the eight curves depicted in Figure 6A indicate a logarithmic relationship 
between the tr«,smission percantage through the optical medhim. and the product of tha absorption coefficient and 
the com^ntration of the «.bsunc. within the medium. As a result of this non«.earity: no representative averege 
spectrum (..e.. one which is substantially path tength Independent) can be obtained. For exampte. , .gnal through 
a path twice as long will typically be attenuated by a factor of 1000 time,, so that'll- an average it, sum would 
35 be inconsequential. 
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If the characteristics of the medium were such that Beer-Lambert's law held firecisety {lb., the curves 
depicted in Figure 6A were perfect exponentials), then a simple logarithm could be used to reduce the transmission 
vector TU) to a set of Gnear equations. Howeverr due to Rayleigh scattering and multiple constituents within the 
finger 130, the intensity of the optical radiation which passes through the fmger are not related according to a strict 
5 application of Beer-Lambert's law. Th» is shown in Figure 6A by a curve 601 which represents an actual 
non-exponential relation between the transmission percentage and the product of absorption, concentration and path 
length. The difference between the curve 601 and the idealized exponential curve 1 is represented by the shaded 
area in Figure 6A. 

In order to transform the relation between transmission percentage and the product of absorption, 
10 concentration and path length depicted in Figure 6 A to a more linear relation, a cubic curve fit is performed on the 
actual transmission vectors to bring them into conformance with the exponential curves depicted in Figure 6A, which 
. is near to taking the logarithm of the transmission percentage. That is, a cubic fit equation is used to bring the 

i 

curve 601, for example, in conformance with the exponential curve 1. Thus, this cubic fit is described hereinafter 
as a "near-log." Over selected regions, a cubic equation may be defined so as to resemtiie, or closely approximate, 

15 a logarithmic curve. Since the actual relation among the values in the transmission vector TU) is typically near 
logarithmic, this near-log (i.e., cubic) equation may be used in place of a logarithm to more accurately linearize the 
values within the transmission vector liA), Of course. It will be understood that a fourth- or higher degree equation 
may be used to provide more accurate linearization of the values within the transmission vector TU). In one 
advantageous ennbodiment of the present invention, the near-log equation has been experimentally determined to be 

20 for transmission through an adult human finger. 

OD - ZBOBB-T^ + 4.7801 'T^ - 3.4215'T + 1.1289 

and is typically applied over the range from T-O to T-1. The precise values for the coefficients in the near-log 
25 equation above may vary from application to application, and, advantageously, an iempirical evaluation with a number 
of test samples is used to precisely define an equation for a given application. 

Once a near-log has been taken of the transmission percentage as discussed above, the relationship between 
the near^log of the transmission and the absorption coefficient times the concentration of the substance withm the 

optical medium is substantially exponential. Thus, by taking a logarithm of the resultant vector values, a Gnear 

r 

30 relationship is observed between the values on each curve as depicted in Figure 68. Although the relationship shown 
HI Figure 6B between the nearlog of the transmission and the absorption coefficient times the concentration is a 
linear one, the slopes of the lines corresponding to different path lengths is more or less steep depending on the path 
length (PI, P2, ... P8). Thus, a different relationship exists between the near-log of the transmission and the product 
of the absorption coefficient and the concentration based upon the path length through the medium. Thus, the path 

35 lengths will not result in a substantially identical output That is, a set of path*length-invariant, i'mear equations has 
not yet been formed. In order to prevent eny error due to accidental differences m path length (i.e., nonlinearity of 
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the values within the transformed TM) matrix), a second log (referred to hereinafter after as the second log) is taken 
of both the corrected transmission percentage and the product of the absorption coefficient and the concentration 
It should be understood here that a near-log fit may also be used in place of the log. as required by the specific 
application to obtain the most foear values. 

The effect of taking the second log is depicted in Figure 6C. As depicted in Figure 6C. each of the lines 
fspresenting diHering path lengths is paraBel so that each of these lines bears substantialy the same rdationsbip 
to the log of the product of the absorption coefficient and the concentration as to the double log of the transmission 
percentage. That is. the values within the resultant linearized optical frequency response matrix OU) are ail Bneariy 
related so that a set of linear equations may be derived from each of the curves depicted in Figure 6C. Finally, aU 
information relating to the pathlength is removed by rotating the curves of Figure 6C so that the curves ere par'aOel 
to the X-axis (not shown). 

The method of Snearizing the optical frequency response matrix is outlined briefly in the flowchart of Figure 
8. As depicted in Figure 8, control passes from a begin block SOO to an activity block 905, wherein the linearization 
process initiates using rned coef fictents that are experimentally detemuned before nin-time. Once enough datd is 
gathered on the patient, to determine that a minimum transmission variance is obtained (due to the pressing motion 
of the finger 130). control of the method passes to an aethrity btock 910 wherein the hear*g cubic curve fit is 
performed on the optical frequency response matrix. That is, eech of the vahies within the transmlssimi vector TM) 
is substituted for T in the cubic equation: 

OD - -2.8088*T' + 4.7801*T' • 3.4215'T + 1.1289 ; 



to form a new. exponential matrix. 

Once a near-hig is taken of the optical frequency response matrix, the first logarithm is taken of the 
transfomied vector values as indicated m an activity block 020. The result of the first logarithm is a matrix having 
25 vahies related to one another as depicted in Figure 6B. 

The second logarithm is taken of the optical frequency response matrix, resulting in the rehition depicted 
m Fflure 6C. as represented in an activity block 930. That is. a logarithm is taken of each of the vah.es within 
the matrix defining the values of Figure 6B to obtain a set of linearly related values. 

A least-squares analysis is performed on the linearized values to detem.ine if the lines defined by the matrix 
(sea Figure 6C) are maxbnally parallel as represented in an activity block 940. That is, the slopes of the lines are 
calculated ami the sum of the squeres of the differences between the cateutated slopes serves as a measure of how 
paraUel the limis are. An Iterative process which involves varying the coefnctonts of the curve fit equation » used 
to determine if the Ones are maximally parayel or at least paraPel enough to obtain a^viible glucose measurement 
Thus, es represented by a decision block 950. a test is performed to determine if the Bnes are parafleL 
If the hnes are not paranel. then control passes to an activity block 955 wherein the vatoes of the coefficients of 
the cubm frt equation are adjusted. In one embodiment, the vahie of one coefficient, begmninfl, for ,x«npte. with 
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the added constant, is increased or decreased by a fixed increment during successive iterations of the adjustment 
routine (Le.r each time the block 855 is entered from the block 950). Control then returns to the activity block 910 
to apply the modified cubic fit equation to the vector values. The process repeats and a new indication is obtained 
of the parallelness of the Bnes. If the lines have become more paralleL then a further adjustment of the same 
5 coefficient is made in the same dvection by the same or smaller increment. This continues until the fines no longer 
become more paralleL If it is determined that the lines become less paraHel based upon an adjustment then the 
same coefficient is adjusted in the opposite d'u^ection (i.e., if the adjusted coefficient was first increasedr it is 
subsequently decreased). This is repeated until a relatively small change is observed in the parallelness of the Enes 
upon successive adjustments of that coefficient. Adjustments are then made to the next coefficient (e.o., the 

10 coefficient on the first order term) in the same manner. Once this process has been repeated for each coefficient 
the decision block 950 determines If the parallelness of the lines defined by the linearized vector is sufficient to 
produce a viable blood glucose measurement 

In. a preferred embodiment an indication is also ghren of the certitude of the glucose determination based 
upon the measured parallelness of the lines of Figure 6C. This measure of certitude may, for example, be displayed 

15 together with the glucose measurement on the display screen 1320 (see Figure 12). If it is determined within the 
decision block 950 that the measurement is not viable, an error signal may be generated; for example, by means of 
a warning light and a new set of measurements is taken to replace the non-viable values within the vector. 

Once the optical density vector has been linearized, the linearized optical frequency response matrix, 0(/l), 
is constructed, as indicated in an acthrity block 960. Control then passes to an end block 970. 

20 OVERALL SIGNAL FLOW 

Figure 9 is a schematic diagram which pictoriatly represents the overall data flow used in accordance with 
the teachings of the present invention to obtain blood glucose concentration by means of optical signal processing. 
As shown in Figure 9, the light source 110 emits light 115 which passes through the lens 117 and the filter 120 
to provide filtered optical radiation 125. The optical radiation 125 passes through the finger 130 to provide an 

25 optical signal 135 used to generate a signal intensity matrix 1000. 

The signal intensity matrix 1000 is multiplied by the Inverse of a filter characteristic matrix 1010 as 
indicated within a block 1005. As shown in Figure 9, the filter characteristic matrix 1010 is derived from an 
analysis of the filter 120, es described above with reference to Figures 4A40 and Figure 7. The inverse transform 
of the filter characteristic matrix 1010 is multiplied by the signal intensity vector 1000 to obtain the optical 

30 frequency response matrix, or transmission vector, 1015. The optical frequency response matrix 1015 is then 
linearized as described ebove with reference to Figure 6A-6B and Figure 8, to obtain the linearized signal strength 
matrix, or linearized optical density vector, 1020. 

The several spectral properties of water, glucose, oxyhemoglobin, deoxyhemoglobin, and scattering as 
represented within the table 1025, are used to construct a blood constituent matrix 1030. The inverse transform 

35 is taken of the blood constituent matrix 1030 to obtain an inverse blood constituent matrix as indicated within the 
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block 1040. The inverse blood constituent matrix is multiplied with the signal strength mau« 1020 to obtain a 
blood concentration times path length vector 1045. 

A ratio is taken between the element of the vector 1 045 indicating glucose concentration times path length 
and the element in the vector 1045 indicating water concentration times path length as indicated within a block 
1050. Thtt mm results in the concentration of glucose in water, which is equivalent to blood fllucose, and is output 
as a data vahie representing glucose concentration. 
METHOD OP PffTOiif^ifir. ppriCAl piypff 

Figure 10 depicts sdiematicail, the general method used in accordance with the present invention to 
manufacturethe optical later 120.1. should be first noted that previous methods employed to fahrl^^^^ 
ftlters typicaUy involved leying out a circular substrate and then selectively increasi«B the coating thlcknes«i. ,„ the 
surface of the circular substrate as the substrate is rotated with uniform speed. 

Such a filter 1500 is depicted in Figure 13 as having coating layers 1502. 1504, 1506, etc.. of increasing 
thmknesses to form a spiral configuration as the fitter 1500 is rotated. Of course, it should be understood that the 
coatan, thickm,sse, depicted in Figure 1 are exaggerated for ease of illustration. This method of optical coating is 
earned around «,b,tafltlally the entire circumference of the circuler substrate so that as the coated substrate 
revolves, the thickness of the optkal coating grows throughout the entire revolution and then suddenly drops beck 
from the thickest coating to the thinnest coating at the end of omi revolution. J 

It has been found, however, that such methods of optical coating require high immn and are extremely 
costly. Furthermore, manufacturing these f ihers is typically carried out omM,y.one. since production methods do mit 
20 eflow for laying out several disks on a single sheet for mass production purposes. 

In accordance with the manufacturing method of the present invention, a flat substrate 1 100 is coated with 
opt«al coating of incraesmg thickness to form wedge-shaped coated layers 1 1 1 0. Of course, it should be noted that 
for purposes of clearly iyustrating the present mvention. the thickness of the optical coating 1110 has been 
eiaggeretad. and in pract»al applications. the thickness of the optkal layer 1110 wiU very from 1.66 micrometers 
to .bout 3.33 micrem«ers, with an overage thickness of about 2.35 micrometers. Of course, it wiO be understood 
that these ttacknesses are approximate and may vary depandinB upon the index of refractton of the layer materials. 

Such a method allows for deposition of a minimal number of optical coating' layers. In one preferred 
embodnnent. only 17 layers are necessary to obtain the desired resolution. In one ertbodiment. alternating tayer, 
of h.gh (2.20) end low (1.48. mdices of refraction are deposited onto the substrate. Although the manufacturing 
method of the present imrention may resuh in less perfect filters than other more Expensive procedures such 
-^perfections can be accommodated in digitel signal processing step, as described above. For example, previous 
filters typ«ally pa«. a single frequency band at a time, white the filter of the preferred embodiment may allow for 
cnultipla bands to pass, since this is accounted for by the signal processing of the invention. 

Once the optical layer 1110 ha. been applied to the substrate 1110. a cylindrical portion 1120 is cut from 
the wedge-shaped sleb formed by the optical layer 1110 together with the .ubstrata 1110. A cylmdrical aperture 
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is then formed in the center of the cylindrical portion 1120 to form the filter 120 and, subsequently, an opticaliy 
opaque strip such as the brass strip 122 is formed over a portion of the optical filter disk 120. 

The above description provides ease of illustration for understanding the essential aspects of the invention. 
However, it should be understood that the method may, in practice, Invohre first cutting the substrata mto a disk 
5 with a shaft; Thereafter, the optical coatmgs are applied onto the disk as though the dislE were still square so that 
the excess faDs onto the platform (not shown) supporting the disk withki the vacuum tank. In this manner the 
wedge is formed on the surface of the disk 120 as shown in Figure 10. 

In one advantageous embodiment, the production specifications for the filter 120 are as follows: 

10 SIZE: 20 mm wide x 20 mm wavelength span, linear multilayer coating 

SUBSTRATE: 25 mm 00 glass disc with 7.5 mm shaft hole in center 

WAVELENGTH: 700-1400 nanometers 

1/2 BW: 50 to 200 nanometers, bands may repeat 

BLOCKING: none 

15 ENVIRONMENT: Survive condensing humidity, 0-70 c 

The pass band edges are produced so as to differentiate a 20 nanometer band edge. 
The pass band may repeat within the window at as little as 400 cm ' spacing, or 17-18 periods within the 
window. The pass band center transmission should approach 100%, and the region between pass bands should 
20 approach 100% reflection. 

Blocking requirements outside of the window are not critical They may be limited by absorption or band- 
edge materials such as RG660, R6700, or semiconductors, or O H bands typically found ;in glass below 7100 cm'\ 
Only the abifity to resolve wavenumber bands near 200 cm*^ with one or more band edges should limit the 
cost. The system is cost sensithre to the filter disc for production quantities of 1,000 to 100,000 per year. 
25 CHARACTERISTICS FOR PRESENT EMBODIMENT 

Preferably, the filter will not have a window narrower than 8,000 to 11,000 cm ^ or about 910 to 1,250 
nm. The bandwidth is advantageously wider than 200 cm \ and the band edge is advantageously narrower than 200 
cm'\ The transmission maximum of the primary band is advantageously above 80%, and the transmission minimum 
is advantageously below 20%. Any other bands should be repeatable, unit to unit; but if they are not, a calibration 
30 ROM could be used in accordance with the DSP to perform initial calibration of mdividuat filters. 
MECHANICAL BOUNDARIES AND CHARACTERISTICS 

The linear filter is advantageously rotated about its center at less than 4,800 RPM, with an aperture 
centered at a radius of minimum 9 mm to maximum 45 mm, with a clear aperture diameter of 1 mm to 3 mm and 
a numerical aperture of .12 to .40. The light path passes through an annular region of the rotating filter, causing 
35 a sinusoidal scan of the wavelengths, although they are deposited finearty. 
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Fof dynamic balance and low turbulence, the linear filter is deposited on a circular substrate. Since the 
center is not used optically, a sundard diameter shaft mounting hole is preferred; most of the present hardware in 
the imrention use either 0.5000-.000. *.00Q5- diameter, or 7.5 -0.0*.1 mm. For a small filter, e.o.. 20 mm 
diameter, bondkiB to the oncoated side would be considered. Note that the filter mount does not have spokes or 
other structural mtemqitian of the optical path. 

Initial optical^anical alignment of the coating on the glass is not critical beyond 0.5 mm and will be 
established electronically. Soma marking of the deposit alignment at the edge or center is desired. 
OVERALL SgnilF WCE OF HAT^ PUn^ '^fffrT 

The foHowing is a summary of the sequence of steps used to process the data signals detected by the 
detector 140. As described above, the sequence of signal processing is extremely important for obtmning an 
accurate meesurement of the patient's blood glucose leweL SignHicant signal processing and analysis is wasted in 
previous systems which did not empkiy the proper sequence. 
THE INSTRUMENT SYSTEM 

The signal generation begins with light source and finger transmission properties. For the glucose sensor, 
start with the receWed transmission, and perform inverse operations on the signal unta only the desired parametei 
is toft Note that the inverse operations are only mrniericaOy equivelent, not the invers^ functhins of mathematicayy 
analytical functions. Each operation is some Uansf ormation which re-orgenizes the, data without convolving or 
otherwise mixing the signal with the noise. The information content of the signal shouU not be decreased or 
corrupted at any stage. Many spectroscopic instruments fail this criterion. Thus, this system design is intended 
to be unique in this field by being in the appropriate order, with a clear objective that avoids corrupting signal 
operations. 

1. TEST TRANSDUCER SIGNALS 

Perform Signal and Noise Classiittation: Synchronous signals correlated in model bounds are deterministic 
(with the same autocorrelation features) and stochastic (random added noise with no aotocorretation features beyond 
25 Gaussian). 

Specificaily, each new received spectrum is added to the deta matrix bhick unti the classifier detecU a 
new noise distribution, or any significant change (detected by means of a statistical test) in the correlatton matrix 
of the 15 non-zero wavelengths. When a change occurs, it indicates that the finger V^O has changed to a new 
position, or that a new finger has been inserted, or that a new patient is in the optical Wth. If hlock-matrix algebra 
« beng used, then at thb change, the system stops adding to the old block and starts creating a new block, et toast 
15 samples tong, to perform statistical tesu on data set. 

If a recurshre method is used, the recursion time constant (sampte number) Is shortened when the classifier 
detects a sample change unti noise Im^eases to a set statistical fimh. or the block rUes the minimum number 
of samples (typically equal to the number of wavetengths resolved. 15 or 16). 

Each answer displayed by the glucose sensor and computer wiQ be the result of one block (of variahte 
length), or of a length set by any recursion process's time constent. Whenever the patient changes finger position. 
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fingers, or even patients, the system observes the change, and starts a new calculation of parameters, including 
glucose concentration. Thus the patient and doctor can have confidence in the readout by knowing every time the 
properties change, the instrument generates fresh data. At the same time, this feedback to the patients encourages 
the patient to remain steady durmg the necessary analysis tone. For the present signal-to-noise ratios, and for the 
5 amount of time a patient is expected to remain steady, the analysis time should be on the order of 3 to 30 seconds. 
When used as a continuous monitor, longer steble periods will be available, and averages on the order of 30 to 180 
seconds might be typical 

2. CLEAN-UP ELECTRONIC NOISE 

Perform the first linear operations to remove electronic noises: subtract linear interference and average 
10 random noise. 

Deterministic noise may be adapted to, and subtracted due to a persistent, periodic form and phase; 
stochastic noise, such as impulses, start-up transients, or white noise spectra, have no'predictable properties. Model 
signals of a priori information become first approximations, and are loaded into any adaptnre or recursive filters, at 
start-up, and are used to predict signals through any impulsive or transient noise. This is the general form of a iow- 
15 pass fitter, but with the feature of averaging the shape of entire waveforms in a heartbeat vector (e.g., a heartbeat 
plethysmographic (equal to pressure, volume) waveform as a vector with one average beat-to-beat period of systolic, 
dicrotic, and diastolic shapes). The noise in this waveform is preferably averaged to zero. 

At the detector 140, ambient light, powerline noise at 60, 120, etc. Hz. eiectrostatic, and magnetic fields, 
and electromagnetic waves are added to the signal photons coming through the patient finger. These noises are 
20 filtered, subtracted or averaged-to-zero before any non-linear stage is encountered. 

When ail linear relationships are removed, only non-linear relatiensh^s remain. The next operation is a non- 
linear transform. Subsequently, ynear operations are performed agam. The last step Is a foiear output of some 
parameter, so the preceding steps are a sequence of linearizations followed by finear operations to uncover signals 
in interference and noise. 

25 The industry standards use light-source carrier*modulation to get away from 60, 120 Hz (50, 100 outside 

USA), and narrow band filtering around the carrier. Such a frequency filter does not remove out of phase noise. 
Dvect subtraction of a best fit noise model with matched phase removes ai) ambient interference. 

3. PREPARE LINEAR SOLUTION 

Linearization, Log Transforms: The object of the operations is to brthogonalize (make stetlstically 
30 independent; make normal; decorrelate) the desired signal from a successkm of added, multiplied, or distorting noises, 
untH only the desired signal remains. 

If the signal is not linear, then it may have multipte vahies for the inverse. Then, no operator or process 
can discriminate the answers. 

4. LINEAR ALGEBRA 

35 Perform matrix algebra on Linear data: Standard spectroscopic matrix algebra solutions are used after ali 

of the assumptions for linear relationships and for stationary signal and noise distributions are satisfied. 
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5. TEST PHYSIOLOGICAL SIGNALS 
Perfom, error detection: determine amplitude, frequency boundaries, and correlate these with any a priori 

hrterferences to continue signal and noise classification. 

6. FEEDBACK 

Feedback to classification. Dnearization or physiological noises: Subtract any Knear interference and average 
the random noise. 

7. FEEDFORWARD 
Apply clinical algorithms and generate derived parameters. 

8. TEST CUNICAL LIMITS, ALARMS 
Define the finger 130 in the optical path. Test if it is mostly water with brge scattering extinction, (about 

30db lossllO mm is typicaHy observed), check the pulse, correlate with previous records to identify the patient and 
the characteristics of the patient's finger. Display "Not a healthy finger" or equivalent, but do not stop operating 
since calibration tests may be done with in vitro materials. If any unusual correlations or signals are present, display 
increasing number of auxffiary parameters as a test, and for user credSiility. 
IS 9. «0. COMMANDS. CUNICAL DATA. DISPLAY. RECORD 

This portion relates to input keys, input data and output format. Ghicose tdkramxs as weU as other 
parameters and historical records ere input. Other embodiments may also require iiiput of diagnosis, predktions. 
instructions per prescription, etc. Afl of the clinical data required to manage a diabetic diet such as the lessons from 
the Diabetic Complications and Control Tests (DCCT), can be eventually put into the program. This data may be 
20 programmed into a personal computer to minimize the sensor computer. 

Although the preferred embodiment of the present invention has been described and illustrated above, those 
skiUed in the art wiH appreciate that various changes and modifications to the present invemion do not depart from 
the spirit of the imrention. For axampte, the principles ami method of the presem invention could be used to detect 
trace motecutes within the btoodstreem (e.B, for drug testing, etc.). In addition, a singte near-log equation of 
sufficient accuracy may be used to linearize the optical frequency response matrix. Accordingly, the scope of the 
present invention is linlted onhr by the scope of the f olknwing eppended claims. 
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WHAT IS CLAIMED IS : 

1. A system for non*invasively monitoring blood glucose concentration within a patient's bloodstream, 
said system comprising: 

a light source which ennts optical radiation at a plwality of wavelengths; 
5 a receptacle for recehring a fleshy medhim of said patient; 

an optical detector positioned to recehre light from said light source and attenuated by said fleshy 
medium, said optical detector responsive to optical radiation of at least said plurality of wavelengths to 
generate an output signal Indicative of the intensity of said optical radiation; and 

a signal processor coupled to said detector to receive said output signal, said signal processor 
10 responsive to said output signal to isolate portions of said output signal due to optical characteristics of 

said fleshy medium to provide a set of frequency response values, said signal processor having a 
linearization module which linearizes said set of frequency response values and to analyze said linearized 
data to determine the concentration of glucose within said patient's bloodstream. 

2. The system of Claim 1, wherein said light source comprises a plurafity of emitters, each emitter 
15 transmitting Kght at a selected one of said phiraOty of wavelengths. 

3. The system of Claim 1, wherein said Gght source comprises a broadband light source, said system 
further comprising an optical filter which selecthrely passes ones of said plurality of vvavefengths. 

4. The system of Claim 3, wherein said detector comprises a single detector responsive to said ones 
of said plurality of wavelengths to provide and output signal indicative of the sum of the Jntensittes of said ones of 

20 said plurality of wavelengths. 

5. The system of Claim 1, wherein said optical detector comprises a plurality of detectors, each 
detector responshre to at least one of said plurality of wavelengths to generate an output signal indicathre of the 
intensity of said at least one wavelength. 

6. The system of Claim 1, wherein the Gnearization module comprises a double logarithm operation. 
25 7. A system for non invasiveiy monitoring blood glucose concentration within a patient's bloodstream 

via spectroscopic measurement of attenuation of a fleshy medium containing blood, said system comprising: 
a light source which emits optical radiation at a pluraOty of wavelengths; 
an optical filter which setectwely passes ones of said plurality of wavelengths; 
an optical detector positioned to receive fight from said fight source which has been fOtered by 
30 said f ater and attenuated by said fleshy medium, said optical detector responsive to fight of at least said 

ones of said plurality of wavelengths to generate an output signal indicative df the intensity of the optical 
radiation filtered by said filter and attenuated by said fleshy medium; and _ 

a signal processor coupled to said detector to recehre said output signal, said signal processor 
responshfe to said output signal to isolate portions of said output signal due to optical characteristics of 
35 said fleshy medium to analyze said portions to determine the concentration of glucose within said patient's 

bloodstream. 
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8. A non-invasive blood glucose monitoring system which analyzes blood glucose within a fleshy 
medium containing blood comprising: 

a light source; 

an optical detector responsive to light from said fight source to generate an output signal; 
V ' compression device which is configured to provide for physical perturbation of said flashy 

medhim to express flind from said fleshy medium; and 

a signal processor responshra to an output signal from said optical detector when fluid is expressed 
from said flashy medium and responshre to an output signal from said optical detector when fbid is not 
expressed from said fleshy medium to isolate information relating to the concentration of glucose in ««d 
■O blood. 

9. A method of non-invasively determinino blood glucose concantration ciomprising the steps of: 
generating a set of vah<es indicative of optical characteristics of significant blood constituents; 

illuminating a fleshy medium having blood with fight of a pluraOty of wavelengths; 
detecting said Kght after attenuation of said light by said fleshy medium; 

" "•""'"'^ " "Bnal from said detected light, said signal indicative of optical characteristics of smd 

fleshy medium; 

isolating components of said signal which are indicative of the concentration of glucose in water 
within said blood in response to said detected light and said set of values indicathre of optical 
characteristics of said significant blood constituents; and 

20 

generating a value indicative of the glucose concentration wHhin said blood. 

10. The method of Claim 9, wherein said significant blood constituents comprise water, hemoglobin, 
oxyhemoglobin end glucose dissolved in fluid. 

n. A method of non-invashrely determining a concentration of a blood constituent comprising the steps 

of: 

'"^"""8 a medium having blood with optical radiation, wherein the^blood has a concentration 
of the blood constituent and the optical radiation traverses a path length through the medium; 

detecting the optical radiation to produce electrical signals indicathre of optical characteristics of 

the medium; and 

linearizing the electrical signals to provide an indication of the concentration of the blood 
30 constituent. 

12. The method of Claim 1 1. wherein said indication of the concentration does not very significantly 
with respect to the path length through the medium. 

13. A method as defined in Cbiim 11. wherem the blood constituent comprises blood glucose 

14. A method as defined to Claim 11. wherein the inearizmg step compriies the steps of- 
in response to said electricel signals, genereting a first set of vahies imJicative of the optical 

characteristics of the mednmi; 
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taking a logarithm of the first set of values to generate a second set of values; and 
taking a logarithm of the second set of values to obtain a linearized set of values that are 
indicative of the concentratnri of the blood constituent. 

15. The method of Claim 14, wherein said concentration of the blood constituent does not vary 
significantly with respect to the path length through the medium. 

16. A method as defined in Claim 11, wherein said linearizing step further comprises the step of 
transforming the first set of values using a polynomial equation to provide a transformed first set of values. 

17. A method of analyzing a concentration of a blood constituent from an optical signal incident upon 
a fleshy medium having blood, the method comprising the steps of: 

detecting the optical signal after attenuation from said fleshy medium to generate a signal having 
a mdicative of optical characteristics of the fleshy medium; 

linearly subtracting components of the signal quantity due to electrical noise; 

removing components of the s^nai due to non>blood constituents of the fleshy medium; 

Bnearizing the signal to generate a set of linearly related values indicative of the concentration pf 
the blood constituent; and 

sohnng the linearly related values for the blood constituent concentration. 

18. A method as defined in Claim 11, further comprising the steps of: 

filtering the optical signal prior to the detecting step with a filter having known spectral 
characteristics; and 

deconvolving components of the electrical signal due to the filter prior to the step of solving for 
the blood constituent concentration, 

19. A method as defined in Claim 18, wherein said deconvolvmg step further comprises the steps of: 
generating a set of values indicatnre of the spectral characteristics of the filter to obtain a filter 

characteristics matrix; 

taking an inverse transform of the filter characteristics matrix; and 

multiplying the inverse transform of the filter characteristics matrix with.a matrix indicative of the 
optical characteristics of the fleshy medium. 

t 

20. A method as defined in Claim 17, wherein said solving step further comprises the steps of: 
generating a spectral library containing values indicative of spectral characteristics of components 

of the blood including the blood constituent; 

generating a set of linear equations having coefficients based upon said values indicative of the 
spectral characteristics of components of the blood; and 

solving the finear equations to obtam a value indicathre of the concentration of the blood 
constituent. 

21. A method as defined in Claim 20, wherein the blood constituent comprises glucose. 
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22. A method as defined in Claim 20, wherein the step of generating the spectral Rbrary tnchides the 
step of observing the spectral characteristics of a known concentration of glucose dissolved in water or blood to 
determine the spectral characteristics of glucose. 

23. A method of manufacturing a rotating optical filter for use in blood glucometry, said method of 
manufacturing comprising the steps of: 

producmg an optical substrate having a top and a bottom; and 

depositing layers of optical coatings on said top such that said layers vary in thickm^sses over 
said top of said substrate in a first direction and such that said layers remain substantially constant in 
thickness in a second direction perpendicular to said first dvectton. 
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